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ABSTRACT 

We present spectroscopic observations of the quasar PG 1512+370, covering 
the H/3 hne spectral range and collected at moderate resolution (2-7 A FWHM) 
from 1988 to 1996. The observations show that the blue wing of the H/3 broad 
profile component has changed significantly in flux and shape between 1988 and 
1990 and between 1995 and 1996. A displaced blue peak on the H/3 profile, 
visible in 1988, but not in the 1990-1995 spectra, is revealed again in one of the 
spectra obtained in 1996. The blue peak (in both the 1988 and 1996 spectra) is 
centered at Avr ~ — SOOOtsoo km s^^ from the rest frame defined by the narrow 
component of H/3, and the [O ill] AA4959,5007 lines. 

We discuss several conflicting interpretations of the data. We find that the 
variability of the H/? blue wing is consistent with Balmer line emission from 
regions whose motion is predominantly radial, if variations of the blue wing 
are a response to continuum changes. Alternatively, we note that observed 
H/? line profile variations are consistent with a variable line component as in a 
"binary black hole" scenario. More frequent observations of H/3 are needed to 
distinguish among these hypotheses. 



Subject headings: Galaxies: Active 
- Galaxies: Nuclei - Line: Profiles 
Individuals (PG 1512+370) 



- Galaxies: Kinematics & Dynamics 

- Quasars: Emission Lines - Quasars: 
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1. Introduction 

Understanding the structure of the inner parsec of quasars-which includes the line 
and continuum emitting regions-has proven elusive. Reverberation mapping techniques 
have been applied to several low-luminosity radio-quiet AGN. Though they required large 
amounts of telescope time, they set on a firm ground basic aspects of the Broad Line Region 



(BLR), i.e., the stratification of ionization (see [Netzer fc Peterson 19971 Peterson 1993| , and 
[Fenston 1991| for reviews). Other results are unfortunately more ambiguous. One reason 
for the ambiguity resides in our poor understanding of the constraints on the inputs to the 



line formation processes in the BLR ( Baldwin 1997| ). Statistical studies of line profiles have 



in principle the advantage of an approach less dependent on the poorly known physical 
conditions of the BLR. However they are as yet plagued by small number statistics and by 
the difficulty of obtaining a uniform sample, in terms of observational (i.e., signal-to-noise 
ratio, S/N) and intrinsic properties (such as optical luminosity; see Paldwin 1997| for a 
recent review). 

Radio-loud AGN have been associated to an extreme in the distribution of parameters 
of radio-quiet AGN ([Brotherton 1996| [Corbin fc Boroson 1996|) . However, Marziani et al. 



(1996, hereafter M96) have shown that radio-loud AGN with detected superluminal motion, 
and apparent velocity j3app ^ 10 show ClvA1549 profiles with low EW and remarkably 
strong red-ward asymmetry. There is no equivalent phenomenology in radio-quiet objects. 
Also, widely spaced double-peaked profiles appear almost exclusively in radio-loud AGN 
(see [Eracleous 1997| for a critical review, and [Eracleous fc Halpern 199^ ). Statistical 



analysis of Balmer line profiles suggests systemic, kinematic, and structural differences in 
the BLR of (at least some) radio-quiet and radio-loud AGN (M96; [Romano et al. 1996| ). 



There is therefore little evidence that the basic reverberation and profile analysis 
results can be extended straightforwardly to radio-loud or high luminosity AGN. In the 
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handful of radio-loud AGN known to have variable profiles, the profiles are often peculiar 
(i.e., double-peaked), as in the cases of 3C 390.3 ( [Veilleux fc Zheng 1991|) , her twin 3C 382, 
OQ 208 ( pVIarziani et al. 1993|) , and Pictor A (pulentic et al. 1995| ; [Halpern fc Eracleous] 



|1994|) . This is different from the case for radio-quiet AGN, since monitored Seyfert nuclei 
do not show double-peaked line profiles. The most common profile among radio-loud AGN, 
however, is not the extremely broad, double-peaked profile observed in Arp 102B; instead, 
it is regular, single peaked, red- ward asymmetric with red- ward peak displacement (M96; 
class AR,R according to Sulentic 1989). The double-peaked profiles are probably not 



representative of the radio-loud AGN as a class. 

The low redshift (heliocentric z = 0.3707 ± 0.0002; M96) radio quasar PG 1512+370 
(B2 1512+37, 4C 37.43) has been frequently observed at optical wavelengths ( P'Brien~^ 
[Gondhalekar 199"T| [Boroson fc Green 1992| , among others). The H/3 broad profile component 



(H/3bc) of PG 1512+370 belongs to the class AR,R for most epochs of observations, and 
hence it is fairly typical of radio-loud AGN. Variability of broad UV lines ClvA1549 and 
Lya has been reported in 1991 from lUE observations ( P'Brien fc Gondhalekar 1991|) . 
These authors suggested the possibility of moderate H/? line flux variations from published 
data (EW(H/?Bc) appeared to vary by 20 %). 

In this paper, we report the discovery of H/3 line profile variations in PG 1512+370, 
and in particular variation of the H/3bc blue wing and of a secondary "blue peak" from the 
analysis of spectra collected from 1988 to 1996 (§0). We will attempt to explain the H/3bc 
profile and its variations as due to (a) accretion disk emission (§ |3.1| ); (b) radially moving 
gas (§|3.2|), or (c) a binary BLR (§p.3|). 
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2. Data and Results 

We summarize the available data in Table 0. The initial reduction followed the standard 
IRAF procedures, which include bias subtraction, flat fielding, wavelength calibration, 
correction for extinction and flux calibration (also see M96, and references therein for 
further details). The spectra have been corrected for redshift, the individual redshifts being 
determined as the average of the peaks of the Gaussian fits to the upper half of H/3 and 
[O III] AA4959,5007. 

The upper panel of Figure |l| shows the optical spectra for PG 1512+370 in the H/? 
region. The spectra have been continuum subtracted by fitting a low order polynomial 
through wavelength regions uncontaminated by emission lines. Most of the weight in the 
continuum fit has been assigned to: (1) the region around 4200 A as adopted for the quasar 
composite spectrum by Francis et al. (1991), and (2) to the red side of [O ill] A5007. Other 
wavelength regions around 4020 A and, with decreasing weight, around 4620 A have been 
included in the fit when available. The spectra have then been scaled to the same fiux 
of [O III] A4959. We did not choose [O ill] A5007 because of the intrinsically larger errors 
due to the presence of the telluric B band at its red side. The scale factor determined on 
the basis of [O ill] A4959 is affected by the extended appearance of the [O ill] emission 
(diameter 4 arcsecs, Durret et al. 1951 ). On the other hand, the H/3 narrow profile 



component (H/3nc) flux is difficult to measure, because the underlying H/3bc is variable. As 
a result, measured fluxes are thought to be accurate within ±10% at a 3a confidence level. 
The procedure to remove the contribution of FEllopt (which is minimal; see, e.g., [Boroson 



fc Green 1992]) , HeiiA4686, H/3nc and [Olll] AA4959,5007 from the H/5 profile is described 



in M96. Table |] reports measurements taken on the cleaned H/3bc line. 

The only signicant change of total H/3bc flux recorded in our spectra occurred probably 
between 1988 and 1990, when it was ~ 30% of the 1988 fiux. The larger decrease is 
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observed in the blue wing of H/9bc and with a possible smaller decrease in the H/3bc red 
wing. Between 1990 and 1996 the H/9bc Aux remained constant within the observational 
uncertainties. No shape change is noticeable in the 1990, 1992, and 1994 spectra. At these 
times, there is no peak at Avr ^ —3000 km s~^ (the "blue peak") or blue shoulder visible 
and the H/3bc profile can be classified as AR,R according to Sulentic (1989). However, the 
line profile showed remarkable changes between 1988 and 1990 as well as between 1995 and 
1996. 

1. The blue peak is visible in the 1988 spectrum but is apparently absent in 1990, when 
the total H/3bc Aux is also significantly lower. 

2. The blue wing appears more prominent in the 1995 spectrum but depressed in the 
1996 spectrum obtained at SPM, with a blue wing fractional fiux decrease of ~ 20% . 
This variation resembles the one that occurred between 1988 and 1990, save that the 
total H/3bc fiux did not change appreciably in the 1995 and 1996 observations. 

The lower panel of Figure |1| shows the difference between the 1995 and 1996 spectra in the 
H/3 region. The slight depression in the 4925-5025 A region is due to the extended nature 
of [Olll], as different contribution from extra-nuclear [Olll] probably fell within the slit in 
the 1995 and 1996 spectra. To compute the probability that the observed change in line 
profile shape is due to random noise fiuctuation, we first estimated the 3a noise level from 
the observed maximum noise fiuctuations over the average signal in the rest frame range 
5060-5100 A (much larger than the noise correlation length Ac, Noise ~ 4 pixels). We rebinned 
the spectrum (over the wavelength range where we observe a change) with a step ~ Ac,Noise- 
We then multiplied the probability that the change in amplitude for each bin is due to noise 
(which is assumed constant over the H/?bc line wing). We find the probability that the 
blue wing change between 1995-1996 is due to random noise fiuctuations is negligibly small. 
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Similar considerations apply to the 1990-1988 H/3bc line profile difference (not shown in 
Figure |l]) . 

The remaining question at this point is whether the blue peak may or may not be 
variable in radial velocity. The ability to resolve the "blue shoulder" into a "blue peak" 
depends on the spectral resolution of the data. The 1988 and 1996 (Asiago) spectra 
have a resolution of 3 A and 2 A FWHM respectively. The blue peak is well visible at 
Av, ^ -2960t^gg km s"^ (in 1988) and Av^ ^ -3400tl^[} km s"^ (in 1996) with respect to 
the radial velocity of H/^nc- The 1995-1996 difference has a broad Gaussian shape, with 
maximum at A ~ 4814 A. Radial velocity values of the blue peak measured in 1988, 1995, 
and 1996 are therefore consistent. The possibility that the blue peak is present in the 1990, 
1992, and 1994 spectra but not visible because of the lower resolution cannot be ruled out. 
After all, the blue flat top observed in the 1996 San Pedro Martir spectrum is resolved as 
a blue peak in the 1996 Asiago spectrum, collected 30 days earlier, when the blue wing 
appeared depressed. Two interpretations seem possible. 



1. The blue peak is a component of H/3bc more variable than the rest of H/5bc; at 
approximately constant Av^.. In this case, the blue peak (i) may be due to a flux 
redistribution within H/3bc unrelated to continuum changes or (ii) may be a structure 
produced by response to continuum variation. Both the 1988 and 1995 spectra (when 
the blue wing was more prominent) show the highest continuum fluxes. This suggests 
that an enhancement of the blue wing is related to an increase of the continuum 
level. However, as the total H/^bc flux did not change significantly in all our spectra 
obtained between 1990 and 1996, it is also possible, albeit less likely, that we are 
observing a line flux redistribution unrelated to continuum variations (as in NGC 
5548; [Wanders fc Peterson 1996]) . Model implications are discussed in §|3T 



2. If the disappearance of the blue peak in 1990-1994 is genuine, the radial velocity of 
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the peak in the blue wing may have changed to ~ OI^ooq km s ^ in 1990, 1992, and 
1994. In § pl3| we discuss a possible model. 

3. Discussion 

3.1. Accretion Disk Line Emission 

The observed H/3bc profile as a whole is not consistent with the predictions of disk 
models at any epoch of observation. Nevertheless, if we wanted to interpret the secondary 
peak Avr ^ —3000 km s~^ as the Doppler-boosted peak of a disk contribution, we would 
need a disk with rin ^ 600, rout ~ 1000, i ^ 21° following Chen & Halpern (1989). 
Non-conventional disk models (elliptical disks or warped disks; [Eracleous et al. 1995| ; 
Marziani, Ualvani, fc Sulentic 1997]) may be more successful in fitting the H/? line profile 
at a given epoch. However, the variation observed between 1988 and 1990 and between 
1995 and 1996 rules out the possibility of "switched-on" disk emission, since the H/3bc 
change appears to be restricted to the blue wing in both cases (see Figure |l]). The observed 
variation would require implausible changes of the disk parameters from one epoch to 
another. 

3.2. A Radially Moving H/3bc Line Component 

The moderate change in the H/3bc blue wing flux, ^ 20%, without a simultaneous 
response in the red wing, points toward line emission from gas whose motion is 
predominantly radial. Given the observed blueshift and assuming optically thin clouds, 
as it is generally believed to be the case for most H/3 emitting clouds, the gas must be 
moving outwards. We are not saying here that the H/9bc line emission is due exclusively 
to outflowing gas, but that a fraction related to the blue peak most probably is. Without 
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further data, it is however not possible to elaborate on this result. Caution is needed as 
the temporal sampling of our spectra is quite probably not sufficient to constrain line and 
continuum variations ( Peterson 19931) . The H/?bc line luminosity is L(H/3bc)~ 1-2 x 10^^ 
ergs s~^ for Hq = 75 km s~^ Mpc~^, and qo = 0.1. Under the assumptions of Menzel-Baker 
recombination theory, the light travel time across the BLR is tlt ~ 107 L^^IJif])^ n^iff^^ 
days. Spectra obtained at a frequency of ^ 1 month^^ are desirable. If the time scale 
for continuum variation is found to be comparable to tlt, then "ripples" in the line profile 
could be produced by reverberation effects. 



At 20 cm, on a FIRST map ( [Becker, White, fc Helfand 1994| ), PG 1512+370 appears 



as a double-lobed radio source, with lobes extending up to a projected linear distance of 
~ 400 kpc from the central core, at P. A. ^ 105°. Integral field spectroscopy revealed two 
[Olll] AA4959,5007 emitting regions approximately aligned along the East-West direction 
with the nucleus ( |Durret et al. 199^ ). Our 1996 spectrum also revealed the eastern emitting 



region as a "blob", at Avr ~ —300 km s~^, in excellent agreement with the results of 
Durret et al. (1994). Extended [Olll] AA4959,5007 emitting regions in radio quasars show 
a statistical tendency to align with the radio-source axis in radio-loud, steep spectrum 
quasars ( [Baum fc Heckman 1989 ). Extended [Olll] AA4959,5007 provides some evidence in 



favor of outflowing gas related to the presence of radio ejecta (see the discussion by [Durret 



et al. 199^ for PG 1512+370). We consider it unlikely that the blue peak is due to some 



out-flowing gas that is episodically related to the radio jet. The dynamical time-scale for 
the BLR is probably shorter than 8 years; a blob moving at Vr ^ 3000 km s~^ could not 
have maintained the same radial velocity over that time period. 
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3.3. Binary BLR 

Within the errors in our spectra, the H/3bc profile variations can be ascribed to a 
stationary profile plus a superimposed component of variable radial velocity, contributing 



up to ^ 10% of the total H/5bc line flux. In a binary BLR framework (see, e.g., [Gaskell 



1983] , [Gaskell 1996a| , [Gaskell 1996b|) , we assume that the secondary component is observed 



at Vr ~ km s^^ (in any case at v^ ^ 1000 km s~^) between 1990 and 1994 (i.e., not 
clearly observed because it underlies the H/3nc)- This assumption is justified since the 
H/3bc underlying the H/3nc ^lay be variable (Figure |ip, and a slightly redshifted "shoulder" 
partly underlying H/3nc may be identified in the 1996 (SPM) spectrum. The "shoulder" 
close to H/3nc is less prominent in 1995, when the blue wing was in turn much stronger. 
This interpretation is not fully satisfactory because the 1988 and 1996 Asiago spectra are of 
higher dispersion than the remaining spectra. The blue peak is resolved in the 1996 Asiago 
spectrum but not in the 1996 San Pedro Martir spectrum, because of the lower resolution 
(the two spectra were obtained 30 days apart). The same may be true for the 1990, 1992, 
and 1994 spectra. 

If we accept the blue peak as a component of variable radial velocity, then the 
immediate implication for the orbit of a binary black hole system is that it is highly 
eccentric, and that the binary must have mass ratio very different from unity. From the 
observed-frame core-to-lobe flux ratio R = 0.091 measured on the FIRST map (1.4 GHz) 
we can estimate the angle between the line of sight and the radio axis. It is known that the 
relationship between R and orientation is valid in a statistical sense, probably because of 
a large intrinsic scatter in the the Lorentz factor of the jet (2 7 20; [Padovani fc Urry 



1992[ ; [Ghisellini et al. 1993[ ), and in the ratio f between the intrinsic jet luminosity and the 



unbeamed luminosity ( Padovani fc Urry 1992[) . We can write R as R = f 5^, where 6 is the 
Doppler factor, and p = 2.75 for a continuous jet ( phisellini et al. 1993| ). We assume 7 =5, 



- 11 - 



and we use f = 7 x 10 ^ ([Padovani fc Urry 1992 ), as appropriate for Fanaroff- Riley type II 



radio sources. We obtain an inclination i ^ 21°. For more extreme values of 7, namely for 
7 = 10 and 7 = 2, we obtain i=15° (f= 7 x 10"^) and i= 26° (f=l x 10~^) respectively. 
To assume i ~ 21° seems therefore reasonable. PG 1512+370 may have been excluded 
from superluminal motion searches, since its radio emission is lobe-dominated. However, at 
i ^ 21° and 7 ?5i 5, we expect that PG 1512+370 should be superluminal with apparent 
transverse velocity /?app ~ 4 (also if 7 = 2, which implies i=26°, we should observe /?app ~ 
1.7). 

If i = 21°, an acceptable fit to the six observed radial velocities {xl ~1.76) is obtained 
for binary period P = 8.15 yr, eccentricity e ~ 0.75, and semi-major axis a = 1.6 x 10^^ cm. 
This solution corresponds to a total mass for the binary of ~ 1.7 x IO'^Mq. The evolution 
of massive black hole binaries orbital parameters (semi-major axis, eccentricity) for several 
mass ratios has been studied recently by [Quinlan (1996)| . From Figure 8 of his work we 



infer that e ~ 0.75 for a = 1.6 x 10^^ cm is somewhat too large for any value of the initial 
eccentricity of the binary, save possibly for the case of initial eccentricity e = 0.9. In this 
case, PG 1512+370 should be in a very-short lived phase of the binary evolution: high 
eccentricity drives large gravitational radiation losses that should lead to the merging of 
the two black holes in ~ 10^ yr. For the lower power Fanaroff- Riley type I radio source OJ 
287 the case for a double black hole is stronger (e.g. |Sillanpaa et al. 1996| ; [Benitez et al. 



1996|) , and rather similar orbital parameters are found: e = 0.678 ± 0.004, P = 12.07 ± 0.01 



yr, i = 4° ± 2° ( |Lehto fc Valtonen 1996| ). If the blazar nature of OJ 287 and the absence 



of obvious line emission is due to Doppler-boosted continuum emission, then PG 1512+370 
and OJ 287 may both host a binary black hole whose orbital plane is observed at different 
inclination. 
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4. Conclusion and Open Issues 

We have presented evidence of moderate variations ( ^ 30% total line flux, with 
the largest variations occurring in the blue wing) in the H/3bc of PG 1512+370 from the 
analysis of spectra obtained from 1988 to 1996. A peak displaced by Avj. ^ —3000 km s~^ 
with respect to the radial velocity of H/3nc (the "blue peak") has been observed in 1988 
and 1996. We have attempted to interpret the observations according to a number of 
presently fashionable BLR models: accretion disk emission (§ ^.11 ); radially moving gas, 
and, in particular, emission due to gas trapped by the radio jet (§pj.2D; a binary BLR 
(§ |3.3| ). Only radially moving gas (but not occasionally associated to the radio jet) and a 
binary BLR are probably appropriate. This is one of the few cases where evidence from 



line proflles may support the presence of a binary BLR (with the cautionary note of § |3.3| ), 
although the eccentricity inferred from the observed radial velocity curve may be too high 
for the estimated major axis. Like OJ 287, PG 1512+370 would then be a very peculiar 
object observed in a very short-lived phase of its existence. However, we favor the idea 
that if H/5bc variations are the response to ionizing continuum changes, the presence of the 
blue peak and the observed variations of the H/3bc blue wing may be indicative of radial 
motion within the BLR. The blue peak may hence be a "ripple" produced by light echoes. 
Indeed, we expect that the blue peak will be detected by future observations at the same 
Avj. ^ —3000 km s~^ therefore ruling out the binary black hole model altogether (as in the 



case of 3C 390.3; pHracleous 1997| ; [Eracleous et al. 1997|) . A set of high S/N ratio spectra 



with a good temporal coverage is necessary in order to test the proposed hypotheses. 
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4600 4800 5000 5200 



A(A) 

Fig. 1.— Upper Panel: optical spectra of PG 1512+370 obtained during 1988, 1990, 1992, 
1994, 1995 and 1996 (see Table |I| for the observation log). Specific flux is in units of 
6.9 X 10~^^ erg s~^ cm~^ (i.e., the average flux of [Olll]A4959 used as a scaling factor). 
The spectra have been vertically offset for clarity, with the exception of the spectrum labeled 
"1996" . Lower Panel: difference between the 1995 and 1996 spectra. Units are the same as 
for the upper panel. 
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Table 1. Data Sample for PG 1512+370 



Observatory Dates Tel. Ap. Spectr. Exp. Time Slit Width P. A. Sp. Res. Detector Ref. 







(m) 




(s) 


(arcsec) 


(°) 


(A) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


La Palma 


1988 May 5 


2.5 


IDS 


8000 


1.6 


a 


3 


235 camera 


A 


KPNO 


1990 Sep 20 


2.1 


Gold 


2700 


1.5 





6.5-7 


TI 800 X 800 


B 


KPNO 


1992 Jul 8 


4.0 


... 


1200 


1.7 


90 


6 


Tek 1024 x 1024 


C 


Calax Alto 


1994 Apr 11 


2.2 


B&Ch 


4800 


2.0 


90 


5.6-8 


Tek 1024 x 1024 


D 


Calar Alto 


1995 May 1 


2.2 


B&Ch 


5400 


1.6 


90 


6 


Tek 1024 x 1024 


E 


Asiago 


1996 Apr 20 


1.8 


B&Ch 


7200 


2.0 


90 


2 


Thomson 7882 


E 


San Pedro Martir 


1996 May 20 


2.2 


B&Ch 


8400 


2.1-2.6 


90 


6-10 


CCDTEK 


E 



Note. — °; not recorded. 

References. — A: Jackson, Penston & Perez 1991; B: Boroson & Green 1992; C: Eracleous & Halpern 1994; D: M96; E: this work. 
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Table 2. PG 1512+370 Rp Line Measurements 



JD EW (H/3) R(3 Flux" H/3 Flux" H/3 Flux" Cont. Specific Flux^ 





(A) 




Blue Wing 


Red Wing 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


2447287 


94 


6.72 


3.49 


3.21 


0.73 


2448155 


124 


5.27 


2.40 


2.87 


0.43 


2448812 


105 


5.23 


2.44 


2.79 


0.49 


2449454 


87 


4.76 


2.52 


2.24 


0.58 


2449839 


65 


5.18 


2.87 


2.31 


0.71 


2450194*= 


75*= 


4.65*= 






0.49 


2450224 


120 


5.03 


2.37 


2.66 


0.53 



Note. — ": H/3 flux in units of [OIII]A4959 flux; ^: continuum speciflc flux in units 
of 10^^ ergs s~^ cm~^ A.~^; the H/3 profile in the Asiago spectrum is truncated at 
4750 A. The EW and flux values reported are lower limits. 



